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Abstract 

Recently, it has been shown that a four-dimensional (4_D) Brans-Dicke (BD) theory with an 
effective matter field and a self interacting potential can be achieved from the vacuum 5D BD 
field equations, where we refer to as a modified Brans-Dicke theory (MBDT). We investigate a 
generalized Bianchi type I anisotropic cosmology in 5D BD theory, and by employing the obtained 
formalism, we derive the induced-matter on any 4D hypersurface in the context of the MBDT. We 
illustrate that if the usual spatial scale factors are functions of the time while the scale factor of 
extra dimension is constant, and the scalar field depends on the time and the fifth coordinate, then 
in general, one will encounter inconsistencies in the field equations. Then, we assume the scale 
factors and the scalar field depend on the time and the extra coordinate as separated variables in 
the power-law forms. Hence, we find a few classes of solutions in 5D space-time through which, 
we probe the one which leads to a generalized Kasner relations among the Kasner parameters. 
The induced scalar potential is found to be in the power-law or in the logarithmic form, however, 
for a constant scaler field and even when the scalar field only depends on the fifth coordinate, it 
vanishes. The conservation law is indeed valid in this MBDT approach for the derived induced 
energy-momentum tensor (EMT). We proceed our investigations for a few cosmological quantities, 
where for simplicity we assume the metric and the scalar field are functions of the time. Hence, 
the EMT satisfies the barotropic equation of state, and the model indicates that constant mean 
Hubble parameter is not allowed. Thus, by appealing to the variation of Hubble parameter, we 
assume a fixed deceleration parameter, and set the evolution of the quantities with respect to 
the fixed deceleration, the BD coupling and the state parameters. The WEC allows a shrinking 
extra dimension for a decelerating expanding universe that, in the constant scalar field, evolves 
the same as the flat FRW space-time in GR. The quantities for the stiff fluid and the radiation 
dominated universe indicate an expanding universe commenced with a big bang. There is a horizon 
for each of the fluids, and the rate of expansion slows down by the time. The allowed ranges of the 
deceleration and the BD coupling parameters have been obtained, and the model gives an empty 
universe when the time goes to inflnity. 

PACS number: 04.20. - q ; 04.50. - h ; 04.90. + e ; 98.80. Jfc 
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1 Introduction 

The original version of the Kaluza-Klein theory [1, 2] assumes, as a postulate, that the fifth 
dimension is compact. However, a particular variety of the Kaluza-Klein theory which is non- 
compact has also been proposed by Wesson [3]- [6], known as the induced-matter theory (IMT) and 
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also as the space-time-mattcr theory. This theory is based on the Campbeh-Magaard theorem [7]- 
[9], which states that any analytic L>-dimensional Riemannian manifold can be locally embedded in 
a {D + l)-dimensional Ricci-flat space-time. The core of IMT applied to general relativity yields the 
4:D Einstein equation as a subset of a 5D Einstein vacuum equation with an effective or induced EMT 
which contains the classical properties of matter. The route usually taken in the IMT is to adapt 
appropriate coordinates that yield the induced-matter consisted with the observations. A generalized 
approach of the IMT, which relates a vacuum {D + l)-dimensional space-time to a Z)-dimensional 
theory with source, has also been proposed [10]. 

Many exact solutions have been found for 5D Einstein eqiiations in the IMT by allowing the 
metric to depend on an extra coordinate and by regarding effects of this extra dimension as an 
induced- matter in a 4D space-time, see, e.g., Ref. [11]. Some of these solutions have remarkable 
mathematical and physical properties. For example, the solutions which are flat in 5D space-time 
but curved in four dimensions provide suitable descriptions of, ordinary matter in the present and the 
early universe [12, 13], and of waves in the inflationary phase [14]. The new solutions in 5D space- 
time have significant applications and relatively more complicated geometry than 5D Riemannian 
one which may be needed to explain the gauge properties of elementary particles. 

Recently, another generalization of the IMT has been proposed for the BD theory as a fundamen- 
tal underlying theory [15]- [17], though the effective equations and quantities in AD hypersurfaces 
obtained in Ref. [15] is different from those in Refs. [16, 17]. Further application of the former 
approach has been performed in [18], though in this work we employ the latter approach and also 
pursue a D-dimensional version of it in another work [19]. In Ref. [17], the MBDT equations are 
obtained from 5D BD equations in which the induced-matter can still be viewed as being generated 
by a fifth dimension. The effective induced-matter} T^v^^^'', is composed of two parts. One part, 
which we present by T^l^^ depends on the BD scalar field, (p, and its derivatives with respect to the 
fifth coordinate. Another part, which we denote by Tji^'^\ is induced from the extra dimensional 
part of the metric and corresponds to the related one in the IMT. An induced scalar potential is also 
produced geometrically from the 5D space-time. As the BD theory can be considered as a modified 
version of general relativity, one may also regard the MBDT as a modified version of the IMT. 

We present and discuss the exact solutions of a generalized Bianchi type I model in the MBDT. The 
Bianchi models are cosmological models that have spatial homogeneous sections invariant under the 
action of a three-dimensional Lie group, but not necessarily isotropic ones. The Bianchi type I model 
corresponds to the flat spatial sections as a homogeneous universe. It generalizes the flat Friedmann- 
Robertson-Walker model (the simplest spatial homogeneous one), which allows for anisotropy. 

The Bianchi type I model for the vacuum solutions are called the Kasner solutions due to E. 
Kasner who first found such solutions in 1921 [20]. When the EMT vanishes, the Einstein tensor, due 
to the Einstein equation, has to vanish as well, but it is not zero in the BD theory. The BD theory for 
the various Bianchi types I-IX models have been considered in the literature, see, e.g., Refs. [21]- [24]. 
Also, the Bianchi type I models with non-vanishing cosmological constant as a function of the scalar 
field, A((/?), have been investigated in the framework of BD theory [25]. In addition, the Kasner 
cosmologies have been considered in the BD theory in four dimensions [26, 27]. Various vacuum 
solutions as well as some dust solutions of the Bianchi type I models with cosmological constant 
have been found in the standard BD theory [22]. Furthermore, the Bianchi type I models have been 
studied in the context of IMT [28, 29] and in the quantum cosmology [30, 31]. 

In the next section, we give a brief review of the MBDT. In Section 3, we investigate a generalized 
Bianchi type I model in 5D space-time, derive the corresponding field equations and meanwhile, point 
to an inconsistency that appears due to the assumption made in the literature. In the rest of this 
section, we solve the equations in a general form under the assumption of separation of variables, 
and obtain five classes of solutions in five dimensions. Among these five classes, we are interested 
in a case (the Case II) which is a generalized Kasner model in 5D space-times. Then, solving 

^In this work, we assume that there is no matter in the 5D space-time. 



the corresponding equations leads to a few relations among the generalized Kasner parameters. In 
Section 4, by applying the effective equations achieved in Section 2, we derive the scalar potential, the 
energy density and pressure of the induced-matter and also discuss the properties of some cosmological 
quantities. Finally, we draw conclusions in Section 5. 



2 Modified BD Theory From 5D Vacuum One 

In the present section, we launch a brief review on the MBDT, the core idea of Refs. [15]- [17] 
while focusing on the latter ones approach. The action of 4D BD theory [32] with a self interacting 
potential^ and a matter field is 



, r(4) 
"T -^matt 
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where^ u > -3/2, c = 1, the D a is the covariant derivative in 4D space— time and the Greek indices 
run from zero to three. The BD scalar field plays a role analogous to the inverse of the Newton 
gravitational constant, the oj is a dimensionless coupling parameter, the Lagrangian -L^^tt depends 
on the metric and other matter fields except ip, and we have chosen c = 1. The variation of action 
(1) with respect to the metric and the scalar field gives the field equations 
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and 
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where S (^L'^lttV^) = -^V^ T^a} ^9°'^ D'^ = D'^D^. In the standard 4D BD theory, one is 
led to the accelerating cosmos in a number of ways. For example, assuming the coupling parameter 
a; as a function of the time [36], or introducing a time-dependent cosmological term [37] and or 
adding a scalar potential to the Lagrangian by means of an ad hoc assumption [38]. However, by 
employing the Wesson idea for the BD theory, it has been shown [16, 17] that such a scalar potential 
can geometrically be induced due to the extra fifth dimension in the MBDT. In this approach, the 
scalar potential has the desired properties and leads to an accelerating universe. In what follows we 
review the approach of Ref. [17] and then, we employ this formalism in the Bianchi type I anisotropic 
model. 

The action for a BD theory in hD space-time, without a self interacting potential, can analogously 
be written as [39] 
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where uj > —4/3, the Latin indices run from zero to four, the Va is the covariant derivative in 5D 
space-time, R^^\ g^^^ and i^^att analogs of the AD quantities. The field equations from action (4) 
obviously are 
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^In the literature, in order to obtain accelerating universes, inclusion of such potentials has been considered in priori 
by hand. 

^In order to have a non-ghost scalar field in the conformally related Einstein frame, i.e. a field with a positive kinetic 
energy term in that frame, the BD coupling parameter must be restricted to a; > — (-D — 1)/{D — 2), where D is the 
dimension of space-time [33]- [35]. 



and 
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where = VaV . The vacuum case is defined as a configuration where there is no matter, except 
the scalar field, in 5D space-time. Thus, by setting T^^ = and consequently T^^^ = g^^T^^ = 0, 



Eqs. (5) and (6) read 
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and 
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Now, following the idea of the IMT approach, one can extract the effective 41? equations. Suppose 
that the 5D line element can be written in local coordinates = {x^, I) as 



dS'' = gAB{x'-')dx'^dx" = g^^{x^)dxi^dx'' + gu{x'^)dl^ , 
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where the / represents the fifth coordinate. This metric is restrictive, but one limits oneself to it for 

reasons of simplicity. It is also favorable to write the fifth component of the metric as (744 (x*-") = 
eh'^{x'-''), where the h{x''') is a scalar field and e = ±1 indicates the sign of the fifth dimension. 
Assuming the 5D space-time is foliated by a family of hypersurfaces S defined by fixed values of the 
fifth coordinate, where the induced metric on each hypersurface, e.g. Eo(Z = /o); is given by 



ds^ = g^^J{x'',lo)dx^'dx'' = gi.udxf'dx''. 



(10) 



The 4D orthogonal unit vector along the extra dimension is written as 
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Dimensional reduction of Eqs. (7) and (8) yields the effective field equations in AD hypersurfaces. 
Actually, reduction of the fundamental field Eq. (7) and the second field equation (8) on the hyper- 
surface So can be written [17] exactly the same as the usual 4Z) BD field Eqs. (2) and (3). However 
now, the term T^^\ that for clarity we indicate it by t'^^^'^\ is the induced EMT of the effective AD 
modified BD theory and has the form 
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The prime denotes partial derivative with respect to the extra coordinate I. The first part of the 
induced EMT comes from the scalar field and its derivatives with respect to the fifth coordinate. The 
second part of T^^ , which is also induced geometrically from the extra dimensional part of the 
metric, is a modified version of the corresponding one in the IMT approach of Ref. [3]. 

Furthermore, in this formalism, these induced Eqs. (2) and (3) can also be derived from an action 
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similar to action (1). Hence, the induced EMT satisfies the conservation law DaT 
the standard BD theory in the Jordan frame. Also, the V{ip) is an induced scalar potential which 
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is defined (or fixed) by the procedure rather than to be introduced priori by hand. That is, on the 
hypersurface Sq one has [17] 
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Obviously, in the particular case uj = —\ and / as a cyclic coordinate, the induced potential, 
without loss of generality, vanishes. Thus, to investigate a general version of 4Z) BD theory with a 
cyclic extra coordinate, one should assume u ^ —1. However, when the metric depends on the extra 
coordinate, one has, in general, a non-vanishing scalar potential. 

In the following, we delineate a generalized Bianchi type I model, as an application of the above 
formalism in the cosmological context. 

3 Generalized Bianchi Type I Metric in bD BD Theory 

Wc apply a generalization of the Bianchi type I anisotropic model for a 5D space-time (with 
spacelike extra dimension) in a comoving reference frame as 



dS'^ = -df + a^{t, l)dx'^ + b'^{t, l)dy^ + c'^{t, l)dz'^ + h^{t, l)df , 



(16) 



where the t is the cosmic time, {x,y,z) are the Cartesian coordinates, and a{t,l), b{t,l), c{t,l) and 
h{t, I) are different cosmological scale factors in each of the four directions. When one assumes that 
there is no matter in 5D space-time and (p = (p{t,l), the field equations (7) and (8), with a little 
manipulation, give 
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where the dot denotes derivative with respect to the cosmic time, Hi, H2, and H^^ are the 
corresponding Hubble parameters in each of the four directions, respectively. 

Meanwhile let us revisit a specific case that has been employed in the literature. When h{t, I) is a 
constant and the other scale factors are functions of the cosmic time only, then the Einstein tensor, 
, will be a function of time as well. Then, if one insists that the scalar field, tp, to remain a 
function of the fifth coordinate in addition to the time, one will encounter inconsistencies when the 
field equations (7) and (8) are solved. For example, these conditions have been used in Ref. [15] when 
the a = b = c. Indeed in this case, for (f{t, I) = (f){t)x{l) one has, from Eq. (23), 

^ + 3i? (*)(/> + |A|(^ = (24) 

and 

§ + |A|X(/) = 0, (25) 

where, without loss of generality, we have chosen the constant h ash = \. The solutions of the above 
equations are 

<^(t) = ^oe(-^^°^^^^^^)* (26) 

and 

x(0 = Xoe^^^^ (27) 

where the Hq is a constant Hubble parameter, the A is a negative separation constant, the ^0 E^nd xo 
are integration constants. However, from Eq. (22) one has 

Go4 = a;^ + — = 0, (28) 

where the last step comes from the fact that the G04 component vanishes for the assumed metric. 
Substituting the separation form of tp with solutions (26) and (27) in Eq. (28) yields 
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The relation (29) for cj 7^ — 1 gives |A| =0. Therefore, x(Z) will be a constant, that is, ip does not 
depend on the fifth coordinate, and thus Eqs. (12) and (15) indicate that T^^ and the induced 
scalar potential are null contrary to the results obtained in the literature [15]. The exception case is 
00 = —1 which is precisely the value predicted when the BD theory is derived as the low energy limit 
of some string theories [40, 41]. 

Also note that, there is a unique solution for when all the metric components (including the extra 
component, contrary to the above condition) depend only on the time while, the scalar field is a 
function of the time and the extra coordinate [17]. 

Based on the usual spatial homogeneity, let us now solve the field equations (17)-(23) using 
separation of variables as 

ip(t, I) = ipotPH">, a{t, I) oc tPH'\ b(t, I) oc t^H'^, c{t, I) oc m'^, h{t, I) = hot^H'^ 

(30) 



where the Hq and ipQ are constants, the p^'s and s^'s are parameters satisfying field equations (17)- 
(23). Actuahy, it turns out to be easier to work with this ansatz. First substituting ansatz (30) into 
Eq. (23) gives 
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This relation is satisfied for the following five conditions 
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Clearly, mixtures of these conditions may also occur. The solution Case I implies is a constant 
which corresponds to general relativity in five dimensions. Such a solution for the spatially flat FRW 
metric has been investigated in Rcf. [17]. In this work, we purpose to investigate the Case II which 
leads to a generalized Kasner space-time in five dimensions. Actually, in what follows we solve its 
corresponding equations and, in the next section by applying the formalism presented in the previous 
section, we investigate the properties of such a solution on each 4D hypersurface. 

Now, substituting ansatz (30) into the rest of the field equations, namely Eqs. (17)-(22), yields 
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L sl + So + S0S4 + SiS2 + S1S3 + S2S3 /(-2*4-2) (37) 

and 

{uj + l)poSo +P2S2+P3S3 = Pii^o + Si + S2 + 53)- (38) 

For later on convenience, we will set /io = 1 in ansatz (30). Also, in order to compare the solutions 
with the corresponding ones in the IMT [28, 29], one must set the BD scalar field to be a constant. 
By assuming p4 ^ 1 and S4 7^ —1, one can match the coefficients of the left and right hand sides 
of relations (33)-(37). Hence, after a little manipulation and considering relation (38), one finds the 
following constrains among the generalized Kasner parameters*^ 
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Constrains (39)-(43) are the generalized Kasner relations in 5D BD theory, though there are only 
five independent relations among the Kasner parameters, and relations (44) and (45) can actually be 
achieved from the others. 

As it will be shown the case po = has particular properties, and hence, we will highlight this 
case in all the following parts of the work. 

• In the special case pQ = = ,sq, the BD scalar field is a constant regardless of any specific value 
of the BD coupling parameter. Thus, constrains (39)-(43) lead to the corresponding relations 
in bD general relativity, Refs. [28, 29], which in turn, give the 5D analogue of the well-known 
Kasner solution. 

• When the scale factors and the scalar field only depend on the time, for the particular case of 
spatial isotropy, by assuming pi = P2 = P3 = P, from (16), (30), (39) and (44), one gets 

dS"^ = -dt^ + AH'^P{dx^ + dy2 + dz^) + Bh^P'^dl'^, (46) 

,i-3p-p4 A 2(3p + p4)(l-P4)-12p^ 
(p = ifot and UJ = , 

(1 - 3p - P4y 



The relations in the second hne of (32), i.e. the Case II, have been rewritten as relations (39) and (41), respectively. 



Q 



where the A and B are constants. These results are the same as the corresponding ones 
investigated for the spatially flat FRW solutions in 5D BD theory [17]. 

In the limit p4 — > (1 — 3p), one obtains again a constant (p and, in general, an infinite uj. In 
this limitf only for p = 1/2, i.e. p4 = —1/2, the field equations are satisfied, namely one has 

dS^ = -dt^ + Ahidx"^ + dy^ + dz^) + Bh-\]f. (47) 

Solution (47) is a unique solution in 5D Ricci flat space-time for metric (16) in the particular 

case where there is not isotropy in all dimensions. 

In the next section, by employing the formalism obtained in the preceding section for the field 
equations, we derive the induced scalar potential, the energy density and pressure of the induced- 
matter, and hence, discuss the properties of the model. 



4 Reduced BD Cosmology in 4D Space— Time 



For metric (16), the induced metric on the hypersurface Sq is just a part of it, and the non- 
vanishing components of induced EMT (12) are 
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where the induced potential l^((/7) should be determined from (15). Also, by interchanging a i — > b 
and a < — > e in relation (49), one obtains {Sir / ip)T'^2^^^'' and (Stt /ip)T^^^'^\ respectively. As, in 
general, the three components j^'p'^'^) (where i = 1,2,3 with no sum) are not equal, one cannot 
consider the induced-matter as a perfect fluid. 

The induced scalar potential is obtained by substituting the power-law relations (30) into defini- 
tion (15) as 



dV 

dip 



2{uj + l)poP4t-^ + Choh-^P% 



-2s4-2 



So 



where 



( = —USo 3so + 2(si + S2 + S3 — S4 — 1) 



+ 2 



S0(S4 - So + 1) + SiS2 + S1S3 + S2S3 



(50) 



(51) 



However, constrains (39)-(45) make C = for the model. Therefore, by integrating (50) while using 
the generalized Kasner relations, one gets 
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on the hypersurface, where we have set the constants of integration equal to zero. Note that, when 
Lu = —1 and or p4 = 0, the effective potential vanishes. 

• For the special case po = 0, the p^ identically vanishes by (50). Therefore, without loss of 
generality, one can set F = in this case. 

In the following subsections we probe the properties of the induced EMT. 



^We disregard the static case p = 0. 
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4.1 Conservation of Induced Energy— Momentum Tensor 



As mentioned, in the MBDT approach, the induced EMT obeys the conservation law as in the 
Einstein theory. Let us check it for the anisotropic induced EMT. 

The conservation law for the induced EMT in the case of Bianchi type I line element reduces to 
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where p = —T 
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,i(EMT) 



(no sum on i) are the energy density and the directional 



pressure of the induced EMT, respectively. 

Now, from relations (30), (41), (45), (48), (49) and (52) one obtains 



87r(po-2) 



Stt 



(po-2)(p4-l)-(a; + l)p§ 
(P4 - 1) - 2{oj + l)ln (^oC*') 



tPo-2 for po 7^ 0, 2 
for po = 2 



(54) 



and 



P4yo'() 
87r(po-2) 

Stt 



{po - 2)pi - (w + 1)pI 



-2 



Pi - 2{u + l)ln ((^oC*') 



for po7^0,2 



for Po = 2. 



(55) 



Obviously, for i ^ j, one has Pi 7^ Pj, that means the induced EMT is not a perfect fluid. By 
substituting the corresponding induced energy density and pressure from relations (54) and (55) 
for the power-law and the logarithmic forms into (53) and using relations (39) and (40), one can 
explicitly show that the conservation law is satisfied. Indeed as mentioned, this result is expected as 
a straightforward consequent of the formalism presented in Section 2. 

• When pq = 0, using the power-law relations (30), equations (48) and (49), and owing to the 
fact that the induced potential vanishes for this case, one easily obtains 



and 



'<goCP4(p4 - 1) 

Stt 

(polo pm 



Stt t2 • 

These results indicate that the conservation law is satisfied for this case as well. 



(56) 
(57) 



4.2 Properties of The Physical Quantities 

In the rest of this work, for simplicityf we assume that the components of metric (16) and the 
scalar field only depend on the cosmic time, namely ansatz (30) reads as 



ip{t) = KtP° , a{t) oc , b{t) cx , c{t) oc t^^ , 
where the k and m are constants. Now, Eqs. (48) and (49) lead to 

Stt o(EMT) ^ _^ , ^i^) 
^ ° ~ h 2ip 



h{t) = mtP\ 



and 



ip zip 



(no sum on i). 



(58) 

(59) 
(60) 



^Nevertheless, we have confined ourselves to a fixed hypersurface. 



in 



Also, from definition (15) and relations (58), one has 



dV 
^ d(p 

Integrating the above relation gives 



2(a; + 1)k2/p"poP4^^^°" 



■2)/po 



Po — 2 



(61) 



(62) 



on the hypersurface, where we have set the constant of integration equal to zero, and have assumed^ 
Po ^ 0, 2. For the case po = 0, one again has V = 0. Substituting relation (62) into relations (59) 
and (60) gives 



_ ^O(EMT) _ K 

P- - 



K / P4 



Stt \pq - 2 



(Po-2)(p4-l)-(a; + l)Po 



and 



K / P4 



(po - 2)pi -{co + l)p'o tP'-\ 



(63) 



(64) 



Stt \po - 2, 

where we have used relations (39) and (40). 

For later on convenience, one usually defines the spatial volume V, the average scale factor a and 
the mean Hubble parameter H as [42]-[45] 



V = abc, 



a = abc = \/v 



and H^Wh, = IZ = ^. 

3^ 3V a 

1=1 



Hence, by relations (39) and (58) one gets 

a(t)(xt'^ with /5 = (Pi +^32 +P3)/3 = (1 -po -P4)/3, 



(65) 



(66) 



where we exclude P = 0, for it gives a static universe. Let us also define the effective pressure on the 
hypersurface as 



1 

3H 



(67) 



i=l 



where this definition yields the conservation law for the isotropic case. Substituting relation (64) into 
definition (67) gives 



PeS = - 



Stt (po-2) 



(po-2)(l-p4) + (u; + l)pg 

(1 - PO - P4:) 



(68) 



where we have used relation (39). Hence, relations (63) and (68) yield a barotropic equation of state 
on the hypersurface as PeS = wp with 



w 



1+P4 



1 - Po - P4 



(69) 



Note that, as po 7^ 2 therefore, w ^ —1. 

Furthermore, in what follows, we present the properties of a few important physical quantities 
such as the expansion scalar 6, the shear scalar and the deceleration parameter q defined in terms 
of the mean Hubble parameter respectively as 



(70) 



'^EYom now on, we do not investigate the logarithmic potential with po = 2, for it leads to some difficulties when the 
weak energy condition (WEC) is applied. For example, in order to satisfy the WEC, one can set p4 = when po = 2, 
which in turn gives no induced EMT on the hypersurface. 
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2 

and 




<'' = k<'" = ilE^f?-3«'l (71) 



These physical quantities are important in the observational cosmology. 
Substituting ansatz (30) into relations (70)-(72) gives 



e = ^-P^-P\ (73) 

a" = ^ [2p4(l - Po - 2p4) - (3a; + A)pl + 2(1 + po)] (74) 

and 

2 + po + P4 

9 = , (75) 

where again we have used relations (39) and (40). Note that, relation (75) indicates that q / —1, i.e. 
the mean Hubble parameter cannot be constant in this model. 

On the other hand, the special law of variation of the Hubble parameter introduced by Bermann [46] 
gives models with constant values of the deceleration parameter. Such a type of assumption has been 
applied in the literature, see, e.g., Refs. [47]-[50]. This type of variation of the Hubble parameter is 
consistent with the observation. By appealing to this law and hence, by considering the deceleration 
parameter to be a constant, i.e. qo, one gets a{t) = {Cit + C2Y^^'^''^^\ where Ci 7^ and C2 are the 
integration constants. With a proper choice of coordinates and constants of integration, the average 
scale factor can be written as^ a{t) oc ^-^^^^ 

comparing with relation (66), means 

/3 = • (76) 
^ g+1 ^ ' 

For an expanding universe, one must set gr + 1 > 0. Thus, any positive value of q corresponds to a 
decelerating universe, while the negative values of g, restricted to — 1 < < 0, yield accelerating ones 
which are consistent with the recent observed data [51]. Now, we can set the most useful characters 
with respect to a fixed q (as an observational quantity), the state parameter w and the BD coupling 
parameter. Indeed, from relations (39), (66), (69) and (76), one gets 

-l + 2g-3^; -{l + q) + ?>w 

Po = — ; and p4 = — . (77) 

q+l q+l 

Substituting pQ and p4 from (77) into relation (63) gives 



P 



(—) 



3w; — (g + 1) 



{l+w){q + lf 



18w^ + 9u;(l - 2q) - 5(1 + 2q) + Aq^ + (1 + 3u; - 2qfu 



t 9+1 . 

(78) 



• In the special case pq = 0, the scalar field is constant and, regardless of any specific value of 
the w, the BD theory reduces to general relativity. Also, we recover the generalized Kasner 
relations of Ref. [28]. In this case, by substituting (57) into (67) and using relation (39), one 
obtains 

Thus, the effective equation of state holds with 

w = , (80 

1-P4 



'For convenience, we will drop the index zero from the constant qo- 
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where (as mentioned before^) ^ 1, and again w ^ —1 (which, in this special case, is also 
consistent with q ^ — 1). 

In order to satisfy the WEC, one must have p4 > 1 or p4 < 0, where the latter choice indicates 
that the extra dimension shrinks by the time. Note that, when — 2 < p4 < 1, the deceleration 
parameter is positive, and when (p4 < —2 or p4 > 1), it is negative. For the radiation dominated 
universe, by setting w = 1/3 one obtains p4 = —1/2 which by substituting it into relations (66) 
and (56) and using relation (39), gives 

a{t) (X t^/^, p{t) oc t'^ for Peff = ^P- (81) 

For the matter dominated universe, we set w = which gives p4 = —1. Then, from (66) and 
(56), one has 

a{t) oc f/^, p{t) oc for P^s = 0. (82) 

In both of the above eras, one has p4 < 0, thus the WEC is satisfied and the extra dimension 
shrinks by the cosmic time. Also for this case, regardless of any specific value of the u, the 
matter and the radiation dominated epoches evolve the same as the corresponding ones of the 
flat FRW space-time in general relativity, as expected. 

In the following subsections, we probe the properties of two interesting cases of the radiation 
dominated era and w = 1 (known as the stiff fluid), for po ^ 0,2. 



4.2.1 Radiation— Dominated Universe 



The main motivation of the IMT in deriving the field equations with matter sources on any 4D 
hypersurface is to interpret the classical properties of matter (viewed as a manifestation of the extra 
dimensions) rather than solely accepting them as a gift. In particular, when the components of 
the metric are independent of the extra coordinate, one obtains the induced EMT to be a traceless 
one [52], in another words, a radiation-like equation of state. Thus, in this subsection, we purpose 
to investigate such an equation of state in the MBDT. 

Substituting w = 1/3 into relations (77) and (78) gives 



2(9 - 1) 
9 + 1 



P4 



'+1 



and 



3R 



eff 



Kq 



87r(g + 1)2 
Thus, relations (73), (74) and (62) give 



2q{co + 2)-q'^{oj + l)-uj 



t 9+1. 



(83) 



(84) 



G 



(9 + l)t' 



-4^2 + 1) + 2g(4a; + 5) - (4a; + 6) 



2iq + 1)2*2 



and 



9 + 1 



2{lO + l)K9-lg(g - if 



(85) 



(86) 



(87) 



(9 + 1)2 

The p > in relation (84) gives — 3/2 < oj < and positive values of q (which is the case for the 
radiation dominated epoch) J° When the deceleration parameter takes positive values, relation (83) 
implies that the extra dimension shrinks as the time increases. The general behavior of the induced 
energy density at an arbitrary fixed time is shown versus q and co in Figure 1. Also, we have plotted 



1 Q 



0.0 



Figure 1: The induced energy density at an arbitrary fixed time as a function of u and q in tlie range of —3/2 < cj < 
and < g < 4 for the radiation dominated universe. 
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Figure 2: The evolution of the induced energy density in the radiation dominated universe with cj = —1.1 and q = 2.5. 
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the induced energy density versus time in Figure 2 for an allowed value of oj and q. 

By assuming the mentioned range that satisfies the WEC for the radiation dominated universe, 
relations (84) (87) show that at t = the physical quantities p, PeS, 6, cr^ and V{(p) take infinite 
values. However, when t — the BD scalar field takes finite values when q > 1 and diverges when 
— 1 < g < 1. As in the limit t — > 0, the spatial volume goes to zero and the expansion scalar 
tends to infinity, thus the evolution of universe starts with zero volume and the rate of the expansion 
is infinite. These results indicate that the model is in accordance with the big bang model. As t 
increases, the expansion scalar decreases but the spatial volume increases, this shows that the rate 
of the expansion slows down by the time. When t — > oo, the spatial volume goes to infinity as well, 
whereas the other physical quantities, i.e. the expansion scalar, the scalar potential, the shear scalar, 
the induced energy density and the effective pressure tend to zero. However in this limit, the scalar 
field has two limits; it goes to infinity when q > 1, while it becomes insignificant when —l<q<l. 
As the ratio fi^/B^ does not depend on the time, the model, in general, does not approach isotropy 
for large values of t. Incidentally, the finite value of the integral 



dt' q + 1 



t 9 + 1 



I , (88) 

to 



to [V(,t')Y/^ Q 
subject to g 7^ 0, shows that there is a horizon in this model. 

4.2.2 Stiff Fluid Distribution 

The high energy density regime of the early universe is usually described with the stiff cosmological 
fluid. This kind of fluid, which is specified by the equation of state Peff = P, has been applied for 
the stellar and cosmological models with utter dense matter [53]. For the stiff fluid, by substituting 
w = 1 into relations (77) and (78), one gains 

2(g - 2) 2 - q 

Po = — Tj-, P4 = ——r (89) 

q+l q+l 

and 

PeS = P= 24!rlq~+l)2 [^^'('^ + 1) " + + (8^^ + H)] (90) 

Substituting these values into relations (73), (74) and (62) gives the physical quantities as 

2 -2g2(a; + l) + (8a; + ll)(9-l) 

" = i^TW' ^ ^ 

and 

2+1 o 

For any value of a; > -1 and < g < (4a; + 7- 3V2a; + 3)/[2(a; + 1)], the WEC is satisfied for the stiff 
fiuid. Figure 3 depicts the variation of the induced energy density at an arbitrary fixed time versus 
the BD coupling and the deceleration parameters. Also, in Figure 4, the behavior of p versus q is 
plotted at an arbitrary fixed time and a constant value of the BD coupling parameter. The evolution 
of p, in this case, is similar to the one depicted in Figure 2. 

For this model, one observes that at f = 0, the spatial volume is zero while the expansion scalar is 
infinite, which shows that the universe starts its evolution with zero volume at the beginning where 



^In obtaining the generalized Kasner relations, we assumed this constrain. 

^°Note that, there are other ranges of w and q that satisfy the WEC, but the mentioned range is the onfy one in 
which the induced energy density behaves smoothty at the whole of the range. 
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Figure 3: The induced energy density at an arbitrary fixed time versus oj and q for tlie stiff fluid, where the BD 
coupling parameter is restricted to — 1 < lj < 100. 
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Figure 4: The induced energy density at an arbitrary fixed time versus q with uj = 0.5 for the stiff fiuid. 
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the rate of the expansion is infinite. At the initial point t = 0, the induced energy density, the effective 
pressure, the shear scalar and the induced scalar potential diverge. However, the BD scalar field has 
no initial singularity for g > 2 and diverges for — 1 < g < 2. By increasing t, the expansion scalar 
decreases while the spatial volume increases, which indicates that the rate of the expansion decreases 
by the time. When t goes to infinity the physical quantities p, P^g, o"^, (p (for —l<q<2), V{ip) and 
O tend to zero, but the spatial volume and the BD scalar field (for q > 2) become infinitely large. 
Thus, in the case of the stiff fluid, one has an empty universe for large values of the time. Again, as 
the ratio cr^/G^ does not depend on the time, the model, in general, docs not approach isotropy for 
large values of the cosmic time. The integral (88) also shows that there is a horizon for the stiff fluid 
distribution, and relation (89) indicates that, for q > 2, its fifth dimension shrinks by the time. 

5 Conclusions 

As the BD theory is a modified version of general relativity, one can view the MBDT as a modified 
version of the IMT in which general relativity is replaced by the BD theory as the fundamental 
underlying theory. We purpose to investigate a generalized Bianchi type I anisotropic cosmology in 
5D BD theory. Then, by employing the obtained formalism, we derive the induced scalar potential, 
the energy density and pressure of the induced-matter on any AD hypersurface in the context of 
the MBDT. Hence, we probe the properties of the induced quantities for the specified cosmological 
model as the quantities fixed by the procedure rather than introduced priori by hand. Indeed, the 
core of the IMT idea, and hence the MBDT, for inducing 5D equations on any 4D hypersurface has 
been motivated to interpret the classical properties of 4D matter rather than solely accepting them. 
Hence, in the cosmological applications, what one obtains from vacuum 5D equations on any 4D 
hypersurface is actually the 5D quantities viewed as the ordinary matter. 

Meanwhile, we illustrate that if a diagonalized metric, in the usual spatial dimensional isotropic 
case and in a gauge with a constant (744, does not depend on the fifth coordinate while the scalar 
field does, it will lead, in general, to inconsistencies when the field equations are solved. Indeed, 
the induced scalar potential and the induced EMT vanish contrary to the results obtained in the 
literature. The exception case is a; = — 1 which is precisely the value predicted when the BD theory 
is derived as the low energy limit of some string theories. Also, there is a unique solution for when all 
the metric components (including the extra component, contrary to the above condition) depend only 
on the cosmic time while the scalar field is a function of the cosmic time and the extra coordinate. 

At first and based on the usual spatial homogeneity, we assume, in general, that the scale factors of 
the specified 51? metric and the BD scalar field depend on the time, and also on the extra coordinate as 
separated variables in the power-law forms. Thus, we find five classes of solutions for the generalized 
Bianchi type I geometry in 5D space-time through which beside general investigations, for the sake of 
compactness, we probe only one interesting class of them that leads to a generalized Kasner cosmology 
in 5D BD theory. This class gives a generalized Kasner relations among the Kasner parameters. One 
important property of these relations occurs when the BD scalar field is a constant. That is, regardless 
of any specific value of the BD coupling parameter, the relations reduce to the corresponding relations 
in 5D general relativity, which in turn, give the 5D analogue of the well-known Kasner solution, as 
expected. Also, when the scale factors and the scalar field only depend on the time, in the particular 
case of spatial isotropy, the results are the same as the corresponding ones investigated before for the 
spatially flat FRW solutions in bD BD theory. 

We then derive the pressure and energy density of the specified induced-matter on any 4D hyper- 
surface. These results indicate that one cannot consider it as a perfect fluid, however we discuss the 
effective AD universe generated by the specified class of the solutions. The induced scalar potential 
is found to be either in the power-law or in the logarithmic form. However, in the particular case of 
constant scaler field and even when the scalar field only depends on the fifth dimension, the induced 
scalar potential, without loss of generality, vanishes. Then, the conservation law (which is supposed to 
be valid in this MBDT approach) has been checked explicitly for the derived induced EMT. Though, 
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as the logarithmic form of the induced EMT has some difficulties in satisfying the WEC, we do not 
investigate it further. 

We proceed our investigations for a few cosmological quantities where for simplicity (and somehow 
without loss of generality) wc assume that the metric and the BD scalar field are only functions of 
the cosmic time. As usual in the Bianchi type I theories, wc define a few (mean) parameters and 
important physical quantities, however, we have defined the effective pressure in a way to yield the 
conservation law in the case of isotropy. The other quantities under consideration are the spatial 
volume, the average scale factor, the mean Hubble parameter, the expansion scalar, the shear scalar 
and the deceleration parameter. First we obtain these quantities in terms of the generalized Kasner 
parameters. Then, we find that the induced EMT satisfies the barotropic equation of state Peff = wp, 
where w; is a function of the Kasner parameters. Also, the model indicates that a constant mean 
Hubble parameter is not allowed. Hence, by appealing to the special law of variation of the Hubble 
parameter, we assume that the deceleration parameter to be constant. And thus, we set the evolution 
of all the quantities with respect to the fixed deceleration parameter q, the state parameter w and 
the BD coupling parameter. In general, the average scale factor indicates that any positive value 
of q corresponds to a decelerating expanding universe, while the negative values of q, restricted to 
— 1 < g < 0, yield accelerating expanding ones. 

In the special case of a constant scalar field, regardless of any specific value of the BD coupling 
parameter, the BD theory reduces to general relativity, and we recover the generalized Kasner rela- 
tions achieved before in the literature. The model, by satisfying the WEC, allows a shrinking extra 
dimension for a decelerating expanding universe which includes the radiation dominated and the mat- 
ter dominated epoches that evolve the same as the corresponding ones of the flat FRW space-time 
in general relativity. 

We then probe the quantities, in general case, versus t and uj for the stiff fluid and the radi- 
ation dominated universe. The results show that in the both fluids, one has an expanding universe 
commenced with a big bang, and there is a horizon for each of them. The rate of expansion slows 
down by the time. Also, by applying the WEC, the allowed (or the well-behaved) ranges of the 
deceleration and the BD coupling parameters have been obtained for each of the fluids. The behavior 
of the quantities, in the very early universe and the very large time, has been discussed. In particular, 
the general behavior of the induced energy density, at an arbitrary fixed time and also its evolution, 
has been depicted in a few figures. The models give empty universes when the cosmic time goes to 
infinity. However, as the ratio does not depend on the cosmic time, the models, in general, 

do not approach isotropy for large values of the cosmic time. 

As presented, the modified 4Z) BD theory with a non-vanishing scalar potential and the induced 
EMT is derived from the bD BD equations with vacuum configuration, contrary to the standard AD 
BD theory where the scalar potential introduced by hand. Note that, though for the Bianchi type I 
model interpreted by the standard approach, introduction of a scalar potential by hand may provide 
more freedom in obtaining interesting solutions, but we have investigated the response of the MBDT 
approach for the specified cosmological model. 
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